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Abstract 

We present a model of a coupled bosonic atom-molecule system, using the re- 
cently developed c-field methods as the basis in our formalism. We derive expres- 
sions for the s-wave scattering length and binding energy within this formalism, and 
by relating these to the corresponding experimental parameters, we can accurately 
determine the phenomenological parameters in our system. 

1 Introduction 

C-field methods have become an indispensable tool in the quantitative description of 
many aspects of Bose-Einstein condensation physics [T|, providing a description of the 
dynamics of highly degenerate Bosonic gases, which incorporates quantum mechanics 
correctly, and is accurate provided the density of the Bose gas is sufficiently high. 

The c-field method is based on the use of a Wigner function representation, in which 
a truncation approximation is made, and this is valid when the density of the Bose- 
Einstein condensate is large — for details the reader is referred to 1 1 1 . The mathematical 
formulation of the c-field method which results from this is superficially very similar to 
that provided by the Gross-Pitaevskii equation. The principal feature in addition to the 
Gross-Pitaevskii equation is provided by the inclusion of a stochastic representation of 
quantum fluctuations in the initial conditions. In the high density limit the method pro- 
vides a treatment in which quantum and thermal phenomena are correctly accounted 
for. The predictions of c-field calculations can be dramatically different from those of the 
Gross-Pitaevskii equation, as was shown in the treatment of colliding condensates (2}{4), 
in which the c-field method was first introduced. 

The quantum fluctuations correspond to half a quantum of noise added to each de- 
gree of freedom; to avoid ultraviolet divergence it is therefore necessary to restrict the 
number of modes used. In practice this is done by means of a projection of the equa- 
tions of motion into a subspace with a maximum momentum, usually called HA. Such a 
cutoff is necessary for two other reasons: 

1. Numerical computations must always be restricted to a finite number of modes. 
This is most commonly provided by the spatial grid, although, as shown in Q] , this 
must be implemented as a projector into the relevant mode subspace if aliasing is 
to be avoided. This issue is not so important for simple Gross-Pitaevskii equation 



simulations, since the occupations of modes at risk of being aliased are usually 
small, but the in the c-field method every mode has at least a half quantum of 
occupation, and aliasing is definitely an issue of concern. 

2. As in the Gross-Pitaevskii equation, the interactions between the condensed parti- 
cles in the c-field formalism are represented by a localized contact potential. This 
pseudopotential method implies the use of a cutoff, as has been argued by Braaten 
and Nieto (5) , who introduced the requirement on the cutoff A and the interparti- 
cle scattering length a s 

Aa s « 1, (1) 
for the validity of a pseudopotential method. 

As implied by our use of the same notation A, the pseudopotential cutoff and that 
required by the numerical simulations are in practice essentially identical. (In prin- 
ciple there is a difference, discussed in the appendix.) In practice the choice of cut- 
off mandated by the numerical algorithms used satisfies this criterion with a satis- 
factory margin of safety; for example in the simulations done by Norrie et al (2HU 
the cutoff used satisfied Aa s a 0.1. 

Apart from the solitary case of hydrogen condensates, in all Bose-Einstein conden- 
sate experiments the scattering length a s used to describe the interactions arises because 
of the existence of a weakly bound state of two atoms. The use of Feshbach resonances 
can give very large scattering lengths, consequently reducing the binding energy of this 
state to a very small value. This can introduce quite slow time scales, which makes it 
wise to investigate the possible influence of molecular dynamics on Bose-Einstein con- 
densation phenomena. 

The c-field method cannot be directly applied to any explicit description of the rel- 
evant molecular physics, since binding can only be described by an attractive poten- 
tial, for which the Gross-Pitaevskii equation and similar c-field equations have no stable 
solutions, and certainly do not produce molecules. In this and subsequent papers, we 
want to combine the ideas of c-fields with those field theory methods which use an ex- 
plicit "molecule field", as originally introduced in [6][7] . These kinds of models are purely 
phenomenological descriptions of the physics, whose parameters must be determined 
to reproduce the correct experimentally measurable quantities. 

The two relevant parameters for the molecular field model are the binding energy of 
the weakly bound state, and the s-wave scattering length, which can be measured. We 
will therefore develop a method for relating them to the phenomenological parameters 
in our formalism. The effective range, we will argue, is not a useful parameter for char- 
acterizing the physics. 

The formalism outlined in this paper is used in Paper II GO, where we implement 
the mean field theory in the Thomas-Fermi approximation, and the Bogoliubov theory 
for this model of a coupled atom-molecule system and use the latter to investigate the 
excitation spectrum for Bragg scattering from a uniform condensate. In Paper III [9) we 
implement this formalism numerically, performing full simulations of Bragg scattering 
from a trapped Bose-Einstein condensate, as in the recent experiment by [10). The re- 
sults from the Bogoliubov calculations and the full simulations both show that the mea- 
sured effects on Bragg scattering as the scattering length increases are well described by 
our coupled atom-molecule formalism. 
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2 Effective Hamiltonian Method for Bose-Einstein Conden- 
sation 



It is normal when describing Bose-Einstein condensation to use an approximate Hamil- 
tonian for the system of ultra-cold atoms that gives a good description of only the long- 
wavelength behaviour of the system, which is all that is relevant for the observable phys- 
ics. The relevant methods are called "pseudopotential methods" or "effective field theory 
methods", and there is a long history associated with the various formulations of these 
methods l5l fTTlfl3l . In all of these formulations, the underlying philosophy is to find an 
approximate description of the the physics of very low energy particles. The definition of 
"low energy" in practice is that the scattering amplitude for such energies does not differ 
significantly from its value at zero energy. 

It is best to formulate the concepts required precisely, and in a form adapted to the 
study of a trapped ultra-cold gas. In such a system all of the relevant physics involves par- 
ticles with a finite small momentum — this means that we can restrict the description to 
particles with momentum less than a fixed cutoff value hA. Consequently, the maximum 
relative momentum is 2fi,A, and for a given centre of mass momentum Q of any pair 
of colliding particles the allowable values of relative momentum, p, satisfy \p + |Q| < A. 
Thus the description of scattering of any pair of particles depends on their centre of mass 
momentum. If the possible values of Q are themselves rather small in comparison with 
HA, this dependence on the centre of mass momentum is not very important, and it is 
normally ignored. 

The appropriate quantum field theory description in the case that the centre of mass 
dependence can be neglected is essentially that of Braaten and Nieto (5), and can be 
characterized as: 

1. The system of atoms is described by a quantum field operator i/>, defined on the 
low energy subspace specified by a momentum space cutoff A 



# C *) = 7^3«( dka k e lkx . (2) 

[2nY'^ Jo 



2. The Hamiltonian for a trapped Bosonic gas of interacting particles is given by 

H = Jdx \^hx)\-^^ + V a {x^{x) + ^fhx)^\x)^{x)^{x)^, (3) 

where V a is the trapping potential. 

3. The quantity U is the inter-atomic scattering strength, and is related to the s-wave 
scattering length a s by 

4nH 2 a s 

£/= £ (4) 

m{\-2Aa s ln) 

For the method to be useful it is necessary that 2Aa s /n « 1, which will mean that the 
predictions of this Hamiltonian are independent of the cutoff, as long as this is not too 
largeQ In this case, at sufficiently low temperatures, the condensate wavefunction T'O, t) 
is accurately described by the Gross-Pitaevskii equation 

_a¥ H 2 V 2x V 4nh 2 a s 7 

lH — = + V(x)y+ ^l 2 ^. (5) 

at 2m m 



1 In fact, there is no reason to believe that the Hamiltonian [3] is valid unless this condtion is satisfied. 
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Experimentally, this is a well-verified equation. 

If 2Aa s ln is not sufficiently small, the simple relationship between the Hamiltonian 
and the Gross-Pitaevskii equation disappears. The derivation of the Gross-Pitaevskii 
equation from the Hamiltonian involves higher order terms in perturbation theory, which 
the choice of a sufficiently small A implicitly sums. 

It is important to emphasize that, although A is often called an "ultraviolet cutoff", 
unlike such cutoffs in quantum electrodynamics, it has a finite value, and this value must 
be small for the pseudopotential Hamiltonian to be valid. Furthermore, in c-field meth- 
ods, all states have at least half a quantum of occupation, so that all momenta up to hA 
participate in calculations, and the corrections discussed in the appendix may be rele- 
vant. 

3 An Effective Field Method for Molecules 

Let us now introduce molecules into the c-field formalism. Since the Gross-Pitaevskii 
equation cannot produce bound states, it is clear that in some sense these molecules 
must be introduced "by hand". 

3.1 The Molecular- Field Hamiltonian 

A method used by several other groups (^[7l [l4lfl6l , is to add an additional field cf> corre- 
sponding to a molecule state, giving the Hamiltonian: 



Here the parameters have the interpretations 

1. Uaa is the background atom interaction strength, leading to the concept of a back- 
ground scattering length at, g = mU ' aa i '4nh 2 . 

2. V a and V m are the external trapping potential for the atoms and the molecules 
respectively. If the magnetic moment of the weakly bound molecule is twice that 
of the atom, then V m -2V a . However, this is not mandatory. 

3. £ is an energy offset term, which allows for a finite binding energy of the molecule. 

4. The coupling parameter g describes strength of coupling of the process by which 
a molecule is formed from two atoms. 

5. The terms with factors U am and U mm correspond to atom-molecule and molecule- 
molecule scattering respectively. Since the atom field is usually much larger than 
the molecule field in the situations we shall consider, these terms are in most cases 
negligible. 

3.1.1 Phenomenology 

The Hamiltonian (6), it must be emphasized, provides only a phenomenological descrip- 
tion of the physics. In particular, the energy offset e can be viewed as a representa- 




^(xl^WfWfd) + Uam^ (X)^ {X)<p{x)f(x) + 




(6) 
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tion of the effect of a Feshbach resonance, and is related to the binding energy of the 
molecule. However, the actual binding energy must be determined by solving the appro- 
priate Schrodinger equation, and it will depend on the other parameters. Furthermore, 
for each value of e, the other parameters, and in particular g, may be different, and in- 
deed, we will find that the values of g and £ required to fit the measured scattering and 
binding properties are strongly interdependent. These dependencies in practice simply 
mean that when setting up a c-field simulation, one must choose the parameters appro- 
priate to the experimental system under investigation. 

In summary, the sole function of the Hamiltonian (|6) is to provide a practical method 
of implementing c-field theory for molecules in a way that is consistent with measured 
properties of the atom-molecule system. The basic theory to which it is an approxi- 
mation is a Hamiltonian involving only atoms interacting through an appropriate inter- 
atomic potential u{x-x'). The molecular field method is only necessary in this paper 
because the "exact" Hamiltonian cannot be represented inside a c-field theory. 

3.1.2 Momentum Cutoffs 

In this case of a coupled atom and molecule system there will be momentum cutoffs for 
both the molecular and the atomic fields, and these can be expressed in terms of projec- 
tors y a and 3> m that project the wavefunctions onto the low energy subspace below the 
cutoff. Because a molecule is formed from two atoms, the interactions will only make 
sense if the molecule cutoff is twice that of the atom. Thus, if the atom cutoff is A, the 
projectors can be defined as 

y a {k) =9(A-|fc|), (7) 

y m m =©(2A-ifcD, (8) 

where is the Heaviside step function. 

In 17181 we have assumed isotropic cutoffs, but this is not necessarily the case for all 
systems. Indeed, as we shall see in Paper III in experimentally realistic systems the cutoff 
can be highly anisotropic. However, for the work in this paper, the exact properties of the 
cutoff are not relevant, and we therefore assume that it is isotropic. 

4 Determination of Parameters 

In order to determine the relationship between this formalism and reality, the parameters 
and the fields have to be related to physically observable quantities. To do this, we need 
to compute 

1. The scattering amplitude, 

2. The binding energy, 

3. The bound state wavefunction. 

4.1 Schrodinger Equation in the 2-Atom : 1 -Molecule Sector 

If \E) is a state in 2-atom : 1-molecule sector, it has a two-component wavefunction, 
defined by the 2-atom amplitude 

y{Xi,X2) = (0|f (Xi)$(X2)\E), (9) 
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and the 1 -molecule amplitude 

(l>ix) = (0\4>(x)\E), (10) 

where |0> is the vacuum state. The normalization of the wavefunctions is given by the 
condition 



J \y{xi,X2)\ 2 dxidx2 + 2 J \(f>[x)\ 2 dx = 1. (11) 

We write the Schrodinger equation for the atom and molecule wavefunctions correspond- 
ing to the Hamiltonian (6) as 

Ey{Xi,X 2 ) = ^ —^{Xi,X2) + S{x 1 -X 2 )(UaaV{Xl,X2) + g(l){X 1 )), (12) 

2m 

ft 2 V 2 g 
E(p{x) = — -^-<p{x) + e<p{x) + ^ty{x,x). (13) 

The equivalent Shrodinger equations for the momentum space wavefunctions y/ and <p 
take the form 

t h 2 \ u r A r A 

\E-—{p\ + pl)\Vf(pi,P2) = — / dq l I dq 2 8(.q l + q 2 -pi-p2)y/(qi,q2) 
\ 2m j {27iy Jo Jo 

+g(p{p\ + p 2 ), (14) 

I 2{iny Jo Jo 



E — e- 



Am 



(15) 



We transform the coordinate system to one of center-of-mass and relative momenta by 
the substitutions P - p\ + p 2 , P' = q\ + q2, k = (pi - P2)I2 and k! - [q\ - q 2 )/2. In the 
centre of mass frame where P = the equations take the form 



E \yr{k) = dk'y(k') + g( p(0), (16) 



h 2 k 2 \ , , U aa rA 

m r (27i) 3 J 

{E-£) (f,(0) = — f-^ [ dk'ifrik'), (17) 
2{2ny Jo 

where we have redefined y/ as the one component momentum space wavefunction for 
the atomic field. 



The projectors 17181 define a preferred frame; in a frame in which P ^ 0, the range of 
the momentum integrals becomes \ k+ jP\ < A. Since the most important interactions in 
a Bose-Einstein condensate do correspond to very small centre of mass momentum (see 
appendix), this is not a great problem. 

4.1.1 The Yamaguchi Equation 

Substituting the molecule function (T7j into the atom equation Q~6), we obtain 
I Em \ C A 

^-p--fc 2 Ji//(fc) = A(£)J dk'y(k'), (18) 



where 



m ( g 2 /2\ 

ME)= m \Uaa+^ ■ (19) 



(2n) 3 h 2 { aa E-e 
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Equation fl8t is similar in form to Yamaguchi's separable potential equation [17], with 
the important difference that our expression has an explicit dependence on the energy 
eigenvalue E in the X factor. However, the method of solution is essentially unaltered. 



4.1.2 Scattering state solution 

A scattering state will be characterized by a positive energy eigenvalue, so we make the 
substitution E — * h 2 K 2 Im in (18), which can now be written as 



(20) 



f{k) = 8(K- fc) + r2 _ fc * + j Q dk'f(k'), 
where we have defined Xk = X{h 2 K 2 /m). Integrating over fc on both sides gives the result 

Small K Approximation The integral on the right hand side can be evaluated exactly, 
and for K « A approximated thus: 



J -a i IA + K\ K 2 
dk = = -4tt A - 2m 2 K+2nK\n a -An A - 2m 2 K + An — . (22) 

o K 2 -k 2 + m \A-Kj A 



r2 

= -AnA-2m 2 K + 2nKln\ 1 -?-^\x -AnA-2m 2 K + An- 

■1T] 

The atom wavefunction for the scattering state can therefore be written as 



i2>(fc) = <5(*:-fc) ^ ~ ~ — \ . (23) 

* {-\IX K -AnA-2m 2 K+AnK 2 IA){K 2 -k 2 +iT]) 

Similarly, the solution for the molecular wavefunction is given by 



^ [ A dk'y{k') 

■Im-E Jo 



g/2 



[h 2 K 2 1 m - e) (- 1 + X K {An A + 2m 2 K - AnK 2 1 A) ) ' 



(24) 



Effective Range and Scattering Length The scattering wavefunction for incident momen- 
tum K and final momentum fc is 

a K {k) = 8{K-k)- Jl 'f? -, (25) 

2n 2 [K 2 -k 2 + ir]) 

where f{k,K) is the scattering amplitude and is thus given by 

-2n 2 X K 112 1 (A + K]]- 1 

f{k,K) = t = = A-iK+-K\n\ .(26) 

l-XK^dkj^— I 2n 2 X K n n [a-kJJ 

The scattering amplitude for small incident momenta can be approximated by the 
effective range expansion, 

1 1 
f(f^) — ~ (271 

J KcotS-iK ~ -l/a s -iK+r K 2 /2' 
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where S is the phase shift, a s is the s-wave scattering length and ro is the effective 
range of the potential. We expand ([26) for small K using the expansions 



In 



1 

A + K 



[2n) 3 h 2 

h 

2K 



Uaa „ 

2e 



+ 7lK' 



4nh 2 



g 2 l2e 



m g[U aa -g z /2e), 



where (29) is valid as long as 
K « A, 



and 1281 is valid as long as 



h 2 K 2 



« 



2U a 



We can thus use equations 1261 and (27) to identify 



a s = 



Anh z 



Uaa „ 
2e 



2 

+ -A 

71 



4 1 (4nh 2 g/2e 



7i A ;r i m (Uaa-g I2e) 



(28) 
(29) 

(30) 
(31) 



(32) 
(33) 



Background Scattering Length Setting g 2 /e — ► in (32) corresponds to the case when 
the effect of the Feshbach resonance is negligible; ie, only the background scatter- 
ing term is nonzero. This means that the background scattering length is 



a bg 



4nh 2 



■ + -A 



mU aa 71 



(34) 



4.1.3 Nature of Effective Range Expansion 



The effective range expansion of the scattering amplitude has the rather limited range of 
validity — the two conditions (30] |3T} set an upper bound for K, and for the expansion to 
be quantitatively valid, K must lie well inside the region defined by the intersection of the 
spaces determined by these bounds. The situation is illustrated in Fig.[T]for parameters 
appropriate to the problems we wish to study, for which available K- space is limited by 
the bound determined by A, that is, by the dashed circle. 

However, the exact solution (26) for the scattering amplitude possesses singularities 
at K = +A, but is analytic for imaginary values of K with modulus very much greater 
than A. In fact, the bound state solution to the scattering problem is given by the pole of 
(26) on the negative imaginary axis and, in cases we consider, is well outside the region 
bounded by the dashed line. This pole cannot be obtained by making the effective range 
expansion (27) for the scattering amplitude. (In fact in Sectior j4~2l we will find it more 
convenient to determine the bound state properties by direct solution of the scattering 
equation, rather than by direct evaluation of the position and residue of the pole, which 
is fully equivalent to determining the pole and its residue). We will therefore interpret 
the available 7<f-space region as the asymmetric region bound by the cutoff A on the real 
axis, but including the bound state on the imaginary axis, indicated by the shaded area 
in Fig.[U 



8 




Figure 1: The grey region shows the region of validity for the scattering amplitude 1261 in the com- 
plex plane. The upper bounds 13011311 for the individual expansions I28II29I are given as the dashed 
and solid circles respectively. The region of validity for the effective range expansion is obviously 
situated within the smaller of these, while the bound state solution (black square) can be found 
far outside this area. The plot here is for a cutoff value of A = 2.2 x 10 6 and a scattering length of 
a = 900flQ» but we get qualitatively similar results for the whole parameter range used in this paper. 



With the model we have chosen, it is not possible to simultaneously fit the binding 
energy, the scattering length and the effective range, since there are only two fitting pa- 
rameters, g and £. The scattering length and binding energy are the relevant measurable 
parameters for the kind of problem we wish to study, whereas the effective range does 
not play any direct role, and therefore we will determine g and £ by fitting the measured 
values of scattering length and binding energy. In this paper, we will not concern our- 
selves any further with the effective range. 



4.2 Bound state solution 



In this atom-molecule model, a physical bound state will be a superposition of a state 
of two i//-field "atoms", and a state of one 0-field "molecule". A bound state solution is 
characterized by a negative energy eigenvalue, E — ► -tl 2 a 2 /m in (18), where a 2 > 0, so 
that K —> (0,ia). Equation IT81 can now be written as 



(35) 



where we have defined A a by setting E — ► -T?a 2 Im in (19). By taking the integral over 7c 
on both sides, we can remove the dependence on the wavefunction, to get 



1 = -4tt 



Jo 



Kk 2 



■dk = -AnX a 



( A 

A - a arctan — 
I a 



(36) 



)o a 2 + k 2 

so that a, and hence the binding energy, are determined by the solution of the equation 



a arctan 



(- 

\ a 



in 



2n 2 h 2 



U„a + ■ 



g 2 l2 



-h 2 a 2 l 



in - 



+ A. 



(37) 
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Threshold of binding At the bound state threshold, where a — «■ 0, we can solve (37) ex- 
actly in terms of e to give 



^threshold — ^ 



U aa +- 



2n 2 h 2 
mA 



(38) 



which is of course identical to the solution of the scattering length equation (32) 
when a s = 0. Substituting from (34), this takes the form 



^threshold 



2 J \47ih 2 > 1 



7T/2A 



(39) 



■ 2Aa,},g /n 

This is nonzero unless g 2 — ► at threshold. 

Weakly bound molecules For small values of a we can expand the left hand side of (37) 
to give 



in a\ \ m I 
A2 ~a' = [2n 2 h 2 { 



m \ Uaa , S 2 I2 
2n 2 h 2 \ -h 2 a 2 1 m- e 



2 A 

so that we get for weakly bound molecules 



+ A, 



e 2 l2 

Uaa+~ 



Anh 2 



2 

+ -A. 

n 



(40) 



(41) 



It is clear that in the case of weakly bound molecules a « 1/ a s , as expected from 
the relationship between the binding energy and the scattering length for suffi- 
ciently small binding energies. 

Tightly bound molecules For large values of a, expanding arctan(A/a) in (37) yields 



r A A J 
, a 3a J 



2n 2 h 2 



U a 



g 2 l2 



h 2 a 2 /m + £ 



+ A, 



(42) 



so that we instead get 



, me n l h l A 6 g £ (2n l A A U aa 

or ~ -r- + — + £ 

h 2 3m I 3 



(43) 



For sufficiently large values of £ we have that a 2 a -melti 2 , which means that the 
detuning e is approximately equal to the binding energy £j, = h 2 a 2 lm for large 
values of the detuning, as would be expected. 



4.2.1 Atomic Fraction of the Bound State 

The bound state solution for the atom wavefunction is 



V(k) = 



N 



r + 



(44) 



where, if the normalization is chosen to be / d 3 k\y/(k)\ 2 = M, the factor N will be given 
by 



m _ r 



dk 



1 



A £2 

An \ — t. 7r-^dk-2n 



o (a 2 + fc 2 ) 2 Jo (a 2 + fc 2 ) 2 



/ 

Jo 



1 (A) A 

— arctan — ^ - 

a \a J A z + a z 



(45) 



10 





lO" 5 



10° 
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Figure 2: The fractions of the atom pairs (dashed line) and elementary molecules (solid line). The 
values of the length parameters are in this case those for 85 Rb; namely l U g a 10~ 7 , a^g = -2.32 x 
10 -8 , and we have chosen A = 2 x 10 6 . 



Similarly the solution for the molecular state can be obtained by substituting (44) into 
fT71 , giving 



0(0) 



r f dk'xfr(k'), (46) 

a 2 lm + E\ Jo 



2[h 2 a 2 /m + E) 

2ngN [ (A\ 
-=— 5 A - aarctan — 



(47) 



AttN I I (A\ \ 2n 2 h 2 \ 

= -_^[/ flfl ^aarctan^-j-Aj-— — j, (48) 

where we have used (37) to get to the last line. 

The ratio of the probabilities for the atom and the molecule states can then be written 



mol l</>(0)| Z (tt^ 

871 



Patm Sdk\y{k)\ 2 I A 

where 



(x) ( ^ A)2 



' (f(a/A)-^/2Aa &g ) 2 ' 



f(a/A) + A 2 /(A 2 + a 2 ) ; 



(49) 



t{z) = zarctan(l/z)-l, (50) 
lu g - (^) 2/3 . (51) 

The result depends on the four parameters which are all lengths, namely l ug , a\,g, A -1 
and a -1 . It is worth noting that while it is mandatory that \K\/A be less than 1, the same 
is not true for \a\IA, which does not necessarily have to be less than 1. On the contrary, 
in most experiments the binding energy is usually such that \a\ » A, for realistic values 
of the cutoff (cf. Section I4T31 . 

The parameter a~ l determines the spatial extent of the atom-pair part of the wave- 
function of the bound state, and this can only be represented on the spatial grid corre- 
sponding to the cutoff A if a I A « 1. The crossover from the predominantly atom-pair 
wavefunction to the elementary molecule wavefunction corresponds to the molecular 
size becoming smaller than the size of the spatial grid. However, {49} shows that, as well 
as A -1 , there are the other lengths a^g and l ug which also come into play. This means 
that the intuitive idea that the crossover happens at a/ A « 1 is significantly modified, 
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I , 1 ! I , 1 

2 4 2 4 



lla 5 [a Q l \ *10" 3 lla s [a- 1 } xl0_3 

Figure 3: Values in SI units of the parameters a) g, and b) e in the Hamiltonian (6), as calculated 
from (32) and (37) for cutoff values A = 2.2 x 10 6 m _1 (hollow black circles) and A = 1 x 10 7 m _1 
(solid red circles). For (a) only values of e corresponding to A = 2.2 x in 6 are plotted, since the 
values of e for the two different values A are almost indistinguishable from each other. In both 
plots the parameters are calculated using data of the binding energy and scattering length taken 
from [18] . Solid lines are to guide the eye. 



and in Fig.[2] which shows the fractions of the atom and the molecule states as functions 
of a/A, it can be seen that the crossover is closer to a/A » 0.01. 

The crossover is thus explicitly cutoff dependent, but this is purely a technical issue. 
There is no fundamental physics in this crossover; it is determined simply by the techni- 
cal need to have a cutoff. This has to be sufficiently large to represent the spatial features 
under investigation, and also has an upper bound determined by the requirement that 
Aa s be much less than one. However, within this range it is arbitrary. 



4.3 Determination of Parameters 

Of the parameters in the Hamiltonian (6) necessary for our modelling, g, £ and U a a, it is 
only the last that can be directly determined experimentally, since the background scat- 
tering length a bg is given by (34). The coupling g and the energy detuning e, however, 
need to be derived from other physically measurable quantities. 

The s-wave scattering length a s is well known for most condensate systems, and 
for 85 Rb the molecular binding energy E b = h 2 a 2 /m has been measured for a range of 
magnetic fields which covers that used to vary the scattering length by means of a Fesh- 
bach resonance. Thus, a range of corresponding experimental values of a s ,a exists, and 
matching these to the expressions (32) and (37) , it is possible to determine the necessary 
sets of values of g and e for a chosen value of the cutoff A in the form 

th 2 a 2 \ {n-2Aa s ){l-2At{alA)a bg ln) 

i 2m ) Aa s {l+t{a/A))-n ' (52) 
1 8nh 4 a 2 \ [a bg (n - 2Aa s ) - na s ){\ - 2At (a/A)a bg /n) 
{ m 2 I 2Aa s (l+ t{alA)-n 
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(a) xi<r 3 (b) 




2 4 1 2 



Va s la l ] *1°~ 3 U a s [a~ l ] xW ~ 3 

Figure 4: a) The ratio E^Ie plotted for cutoff values A = 2.2 x 10 6 m _1 (hollow black circles) and 
A = 1 x 10 7 m _1 (solid red circles). This ratio is seen to be very close to 1, deviating only at scatter- 
ing lengths more than of about 5 x 10 -8 m a lOOOagohr. The deviation at larger scattering lengths 
depends strongly on the choice of cutoff A, and is also affected by mean field effects; b) Plot of 
a vs. aj 1 ; the value of a must approach aj^ for sufficiently small a. However, the data used 
for these plots is measured in a Bose-Einstein condensate, and the four points corresponding to 
aj 1 < 1 x 10 7 m _1 deviate from this law because of mean field effects specific to that situation. In 
both plots the parameters are calculated using data of the binding energy and scattering length 
taken from [18] . 



In Fig.[3]we show the Hamiltonian parameters as functions of the s-wave scattering length 
calculated in this way for two different values of the cutoff A, using data of the binding 
energy and scattering length taken from [18]. The dependence on A is seen to be quite 
weak. (In applications to c-field calculations, the cutoff is related to the size of the simu- 
lation grid and it is important to choose the grid such that a s and £j, depend only weakly 
on the value of A.) 

Figs.[3] and [4] also show that the parameter g is essentially a linear function of a~ l , 
while energy offset e is almost proportional to the binding energy Looking at Fig.|4] it 
can be seen that Zjj, and e are essentially equal except for the very small binding energies 
which occur at large a s . In fact, since the data show that g is unlikely to be zero at the 
threshold of binding, where a~ l — ► 0, it can be seen from 138t that ^threshold ^~ 0> and given 
that for 85 Rb the background scattering length is negative, the threshold value of e is in 
fact positive. 

In panel b) of Fig.|4] it can also be seen that aj 1 and a become very close, but when 
equality approaches at weak binding, the mean field effect of the Bose-Einstein conden- 
sate used in the experimental measurement obscures the equality, which was also noted 
in GU. 



4.4 Validity of the Effective range Expansion 

Using (52) and (53) we can now estimate the validity of the effective range expansion in 
section l4.1.2l In Fig.[5]we plot the ratio of the effective range expansion KcotS a \l a s + 
tqK 2 12 and the exact expression for ^Ccot<5, given by the inverse of the real part of the 
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I 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

K/A 

Figure 5: Ratio of the effective range expansion and the exact expression for KcotS. The dash- 
dotted purple line corresponds to A = 10 7 and the blue solid line corresponds to A = 2.2 x 10 6 . The 
scattering length is set to a s = 900ao and the effective range is given by equation 1331 . 



right hand side of (26). We can easily see that the effective range expansion is only valid 
up to K a 0.6A. 



5 Conclusion 

The Hamiltonian (6) can thus be used to simulate realistic atom-molecule systems, as 
long as the binding energy and s-wave scattering length are known. It is not important to 
match the value of the effective range, which in this kind of model is determined by the 
value of the cutoff, and has little physical significance. The mapping of these physically 
observable parameters to the phenomenological parameters in the Hamiltonian will be 
dependent on the choice of the cutoff, but the results obtained in simulations should be 
essentially independent of the choice of A. 

In Paper II and Paper III, we will apply the model Hamiltonian we have formulated 
here to the application of c-field methods to: 

1. The formulation (for this model) of the simplest mean-field theory, and in particu- 
lar the Thomas-Fermi approximation to the profile of a trapped condensate. 

2. The study of the excitation spectrum of Bose-Einstein condensates and application 
to Bragg scattering of a homogeneous system; 

3. The modelling of Feshbach-resonance-enhanced Bragg scattering from a trapped 
inhomogeneous Bose-Einstein condensate. 

We will demonstrate that the experimental data from (TO) can be reproduced accurately. 
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A Corrections to the Scattering Length and Binding Energy 
for Non-Zero Centre of Mass Momentum 



In Section fOI we assume that the centre of mass momentum of the system is zero, and 
thus we solve the integral in fL8l on the range [0, A] . In this appendix we are investigating 
the corrections that occur when the centre of mass momentum is non-zero. 



A.1 Scattering length 

To find the corrections to the scattering length, we wish to evaluate the integral in t2D 
1 



dq- 



fa Q K 2 -q 2 + lrj 
where the range JIq is now defined by 
|Q±2q| < 2A, 



(54) 



(55) 



for the centre of mass momentum HQ = pi + P2 = P3 + P4- 

We change the integral to polar coordinates, so that the range 31q is equivalent to 



==> |cos0| <X = min|l, 
Thus, we can write 



Q 2 + 4q 2 + 4qQ cos0 < 4A 2 



4A 2 - Q 2 - 4q 2 



4qQ 



(56) 
(57) 



3 S = 2n I — dq \ 

Jo K 2 -q 2 + ir] 'Jo 



Va 2 -Q 2 /4 „2 

= An I — r r dominll 

K 2 -q 2 + irj * 1 



•f 

Jo 



2dcos6 

4A 2 -Q 2 -4q 2 



4qQ 



(58) 



(59) 



Noting that the point at which (4A 2 - Q 2 - 4q 2 ) l4mqQ = 1 is when q = 2A - Q, we can 
write 



q 2 4A 2 -Q 2 -4q 2 



-QI2 K 2 - q 2 + IT] 



4qQ 



+4n 



A-Q/2 



K 2 - q 2 + it] 



dq, 



(A Q A 2 -Q 2 /4-K 2 /2 

4n -"- + In 

2 4 2Q 



K 2 - (A-Q/2) 2 



Q m K 
+4n -A+ — K+— In 

2 2 2 



K 2 -A 2 + Q 2 I4 
K+A-QI2 



K-A + Q/2 



? A 2 -Q 2 /4-^ 2 /2 
-27rA + tiQ - 2m z K + 2n In 



K z - (A-Q/2) 2 



K 2 - A 2 + Q 2 /4 



+2nKln 



K+A-QI2 
K-A + QI2 



(60) 



(61) 



(62) 



15 



In the case when K « A, we can set K —> in I 
amplitude becomes 



and the expression for the scattering 



/ = 



-2n 2 X K 



1 AO lA 2 -Q 2 /4 
- h h In 

2ji z A,k n 2n ti Q 



A-Q/2 



A+Q/2 



giving us the following expression for the scattering length, 



4nh z 



2,-1 

Van „ 
2e 



n 2n n 



1 A 2 -Q 2 /4 



In 



A-Q/2 



A+Q/2 



For Q « 2A, this can be approximated by 



(63) 



(64) 



4nh z 



in 



Uaa „ 

2e 



1 2 Q Q 2 
ti 2n 6nA 



(65) 



Thus, the consequence of a non-zero Q is a reduction in the effective value of A. How- 
ever, this reduction will have very little impact on the value of the scattering length, since 
in practice A must be chosen to reduce any such effects. 



A.2 Binding Energy 

To find the correction to the binding energy, we wish to solve the integral in (36) 



,= [ dq^— 

JK Q a z + t 



3 b = I dq 3 (66) 

where the range 31q is now defined by 

|Q + 2q| < 2A, (67) 

for the centre of mass momentum HQ = pi + P2 = P3 + P4- 

Following the same method as above, we can write the integral as 



r V / X^C?M q 2 4A 2_ Q 2_ 4q 2 .A-Q/2 q 2 

3b = 4tt , 7^- -dq + 4n\ -JJ—dq, (68) 

Ja-q/2 a z + q z 4qQ Jo a z + q z 



4A 2 -Q 2 + 4a 2 fa 2 + A 2 -Q 2 /4 
jiQ - 2nA + ji In 1 



2Q U 2 + (A-Q/2) 2 

(A-Q/2\ 

+4^A-2H-Q-4n-arctan , 



(69) 



(A-Q/2\ jiQ 7 

S3 -4;rarctan — \+4n (A-Q/2) + „ „ (A-Q/2) +..., (70) 

\ a J a z + A z 

where the expansion on the last line is valid as long as Q 2 « 4a 2 + 4A 2 . This is a reason- 
able approximation, since Q < 2A always, and a » A usually holds. The expression for 
the binding energy corresponding to l37l then becomes 



(A-Q/2\ [ m I g 2 /2 \] 1 I Q(A-Q/2)\ 

«arctan(^)s3[_( f / fla+ -^ 7 __)j +(A _ Q/2) ( 1 + (71) 

The assumed condition Q 2 « 4a 2 + 4A 2 means that the final term on the right hand side 
is very little different from A- Q/2, so that the effect of non-zero Q is to replace A in the 
binding energy condition $37\ by A - Q/2 < A. This reduction of A will have very little 
effect on the binding energy, as can be seen in Fig.|U 
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